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Directly observing a zero energy Majorana state in the vortex core of a chiral superconductor 
by tunneling spectroscopy requires energy resolution better than the spacing between core states 
Ao^/eF- We show that nevertheless, its existence can be decisively detected by comparing the 
temperature broadened tunneling conductance of a vortex with that of an antivortex even at tem- 
peratures T ^ Aq^ /cf- 



Driven partially by the dream of building naturally 
error resistant quantum computers, the study of topo- 
logical phases of matter has become an important topic 
of research yy. The simplest class of topological phases 
of matter that could be useful in this respect are the 
chiral Px^Wy BCS paired systems [2]. There are sev- 
eral physical systems where Px + Wy pairing is believed 
to be realized, including the A phase of superfluid ^He 
[3] (^HeA), the exotic superconductor Sr2Ru04[4 , and 
the V = b/2 quantum Hall state [5l|6]. In addition there 
have been recent proposals to realize Px^Wy pairing in 
cold fermion gases [7 . In these (weak) Px^ipy systems, 
certain types of vortices (quasiparticles in the quantum 
Hall context [8]) are believed to carry zero energy Majo- 
rana fermions [H |9] which are the topologically protected 
degrees of freedom. 

In Sr2Ru04 and ^HeA the vortices that carry the Ma- 
jorana fermions are the so-called half-quantum vortices, 
which can be thought of as a vortex in the order pa- 
rameter of one spin species without a vortex in the order 
parameter of the opposite species [10 . (Note that in spin- 
polarized Px-\-iPy systems, including proposed atomic gas 
realizations or the 5/2 state, there is no half-quantum 
vortex and the full quantum vortex carries the Majorana 
fermions.) 

Let us suppose that in one of these systems, the rel- 
evant Majorana-fermion-carrying vortex has been ob- 
served [11]. The next important step would be to design 
an experiment to observe the Majorana fermion in such a 
vortex [12]. In the case of Sr2Ru04, one obvious experi- 
ment would be an energy-resolved tunneling experiment, 
which measures the local density of states (LDOS) [I3] . 
An observation of a localized mode at precisely zero en- 
ergy would be direct evidence of the Majorana mode. For 
cold atoms, an analogous experiment for observing the 
LDOS would be an energy-resolved local particle annihi- 
lation experiment. For the other realizations of Px^Wy 
order it is not as clear how such an experiment would be 
performed [14 . 

In principle such tunneling experiments could pro- 
vide definitive evidence for the Majorana mode. How- 
ever, in practice they may be prohibitively difficult. In 
the vortex, there will exist sub-gap bound states in 
the core known as Caroli-de-Gennes-Matricon (CdGM) 
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FIG. 1: A vortex pair of the Px + ipy superconductor on the 
sphere, described by Eq. [2j Thin white lines represent the 
current flow. Wide arrows represent the pair relative angular 
momentum, is the pairing range. ^ is the coherence length 
which determines the vortex core size. 



states [151 [16]. The spacing between these bound states is 
typically of order Sc = Aq^ /ep where Aq is the gap (pre- 
sumably on order of the critical temperature) and cf is 
the Fermi energy. Since the experimentally observed tun- 
neling spectrum will be smeared by the temperature, this 
tunneling experiment would naively only have a clear sig- 
nature for T < Sc. Unfortunately such low temperatures 
could potentially be unattainable in any of the proposed 
realizations {Sc ~ TjaK in ^HeA, < O.lmK in Sr2Ru04). 
The purpose of this paper is to demonstrate that the tun- 
neling spectrum retains an unambiguous signature of the 
Majorana fermion at much higher temperatures. The sig- 
nature is found by comparing the tunneling conductance 
peaks of a vortex with an antivortex, the direction of vor- 
ticity being defined relative to the angular momentum of 
the chiral order parameter. 

Bogoliubov de-Gennes (BdG) theory. We consider 
a two dimensional uniform Px + ipy superconductor of 
spinless fermions. The BdG excitations are given by [17] 

where T is the kinetic energy operator, and ei? is the 
Fermi energy. 
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FIG. 2: BdG spectrum En(m)i of the vortex pair on the 
sphere, depicting the CdGM core states. The inset shows that 
their double degeneracies are spUt by weak tunneUng between 
the poles. The state nearest zero energy is the Majorana mode 
of both vortex and antivortex. 



We implement the BdG equation on a sphere of radius 
parameterized by the unit vector ft = (O^cj)). The 
spherical geometry has two important advantages: (i) It 
has no boundaries, which strongly affect the low energy 
spectrum, (ii) In the absence of disorder the azimuthal 
angular momentum is conserved, which greatly reduces 
the computational difficulty of the BdG diagonalization. 

The order parameter field on the sphere is taken to be 
of the following form [HI ITS] 
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which defines Aim^Um'- is the pairing amplitude, the 
pairing range is ^p, and Ip is the Fermi angular momen- 
tum, given hy ep = If{If -\- l)/{2mR^). The functions 
a = cos(^/2) and (5 = sin(^/2)e~*'^ are spinor functions. 
Yqim are monopole harmonics [1£ , where q^l^m are half 
odd integers. The order parameter Ap(f2, ft') acquires a 
27r phase when fl encircles which describes Px^Wy 
pairing. |Ap| keeps the particles within the pairing range 

\n - n'\ - ^p. 

The order parameter field F^^iCt) describes the vor- 
ticity of the pair center of mass 11 = (fi + fl')/2. We 
choose Fvv to describe an antivortex on the north pole 
and a vortex on the south pole, with the direction of vor- 
ticity defined relative to the chirality of the Px^Wy order 
parameter, depicted in Fig. [l] For the vortex pair field. 



we use the analytical form (without self consistency) [20] 



E fLYLim, 



(3) 



L=l,3,5,.. 



which defines /l- Ylm are spherical harmonics, and ^ = 
2eF/(7rAo/cF) is Pippard's coherence length. We take 
< i for simplicity. 

The BdG equation is represented as a matrix in terms 
of 3j symbols as 
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Diagonalizing Eq. ([T]) produces a set of energies 
and corresponding eigenvectors u^^^v\^. By azimuthal 
symmetry, m is a good quantum number. The BdG wave- 
functions on the sphere are 

u^m = E<™r (4) 

I 

In Fig. [2] we depict the BdG spectrum of the vortex pair 
as a function of m. The continuum states above the 
gap \En\ > Ao are extended, while the branch that ap- 
proaches zero is the Px^Wy version of the CdGM core 
states. Their number is of order cf/Aq, and their spacing 
is of order 5c [16^. 

As seen in the inset of Fig. [2| each CdGM state is al- 
most doubly degenerate. The splitting represents weak 
tunneling between the north and south pole core states. 
Indeed, we find that the tunnel splittings decrease expo- 
nentially with the radius of the sphere 8 En ^ e~^^^ for 

The lowest positive energy approaches zero as Eq ~ 
Q-R/^, In the large sphere limit, the wavefunctions 
uo{Q) ~ '^o(^) are equally split between the north and 
south poles. The corresponding BdG fermion is con- 
structed out of two well separated Majorana operators 
in the north and south poles. We have also verified that 
Eq is insensitive to the addition of moderate potential 
disorder [2^. This agrees with previous asymptotic cal- 
culations in the plane which have shown that the Ma- 
jorana excitations are 'topologically protected' against 
perturbations [8, 9 . 
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The asymptotic predictions for the wavefunctions 
uo{r) for core- less vortices in the plane [7] are 



uo{x.) 



Jo{kFr)e ^/^^ antivortex, . . 



Ji(/cFr)e*^e- 



vortex. 



which are valid for both r ^ ^ and r <C ^. Numerically, 
we confirmed that these predictions hold even in finite 
core sizes. 

The physical reason behind the difference in Eq. (|6| 
is that the Majorana wavefunctions are sensitive to the 
sum of vorticity and relative angular momentum. Jq is 
obtained only when that sum vanishes, and this is im- 
portant for the experimental signature we discuss below. 

Local Density of States (LDOS). At zero temperature 
the LDOS is defined as [13^ 

r{E,r) = Y,\^n{r)\H{E-En)Mvn{T)\H{E^En)^ (7) 

n 

where r is the distance from the vortex (or antivortex) 
center. 
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FIG. 3: Zero temperature local density of states of Eq. ([7| 
near the vortex and the antivortex centers. The peaks belong 
to the CdGM states. Notice that the zero energy Majorana 
mode is removed from the origin in the vortex, while it is 
maximized at the origin in the antivortex. 



Fig. |3] shows the LDOS near the cores of the vortex 
and the antivortex for displacements r <C ^, and energies 
1^1 < Aq. The Majorana state can be easily discerned as 
the zero energy peak in both vortex and antivortex cores. 

The other CdGM core states also appear as oscillatory 
peaks, with energy spacing 5^ The difference between 
the vortex and antivortex excitations is apparent: In 
Eq. |6]the antivortex Majorana state is peaked at r = 0, 
while the other CdGM states have nodes at r = 0. In 
contrast, the vortex Majorana state is peaked at half a 
Fermi wavelength away from the origin, while two CdGM 
states are peaked at r = 0. Notice that in the vortex, the 
Majorana state has a significantly lower peak than in the 
antivortex. 



In a tunneling spectroscopy experiment (e.g. Ref. [22]), 
the discrete LDOS spectrum is smeared by temperature 
broadening. The tunneling conductance [23] is defined as 



dl 
dV 



E') 



nE',r), (8) 



where f{E) is the Fermi-Dirac distribution at zero chem- 
ical potential and temperature T. 

In the BCS weak coupling regime, kpS, ^ 1, and there- 
fore Sc could be a very small temperature scale. At mod- 
erate temperatures dc < T < Aq^ the peaks of Fig. [3] 
are smeared on the energy axis (but not on the r axis, 
therefore an asymmetry effect can be observed. 
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FIG. 4: Tunneling conductance of Eq. (|8|), in arbitrary units. 
Ai? is the fermi wavelength. The temperature T = 0.15Ao, is 
about 10 times larger than the CdGM level spacing. 



A typical tunneling conductance is depicted in Fig. [4j 
which shows a central peak at r = 0, £^ = 0, with low 
broad ridges dispersing away to larger r, E. We see that 
the central peak of the vortex is twice the height of that 
of the antivortex. 

This effect is a direct consequence of Eq. (|6|. Under 
temperature smearing the two CdGM peaks at r = of 
the vortex, merge into one large central peak. In contrast, 
only a single Majorana state is responsible for the central 
peak of the antivortex. Since the relevant maximas in 
the LDOS are nearly identical, a ratio of 2 is obtained at 
elevated temperatures. 

Our effect requires having spatial resolution in tun- 
neling conductance better than a Fermi wavelength. 
If dl /dV{r^ E) is convoluted with areal resolution of 
(Jr)^ > A|., the ratio between the vortex and antivor- 
tex peak heights rapidly approaches unity as dr > Xp- 
The ratio is weakly temperature dependent in the regime 
5c<T< Ao. 

In real three dimensional samples, zero bias peaks are 
somewhat suppressed by bulk states and surface imper- 
fections, (Ref. [22 reports 15% enhancement above the 
high voltage background). Nevertheless, it is the ratio 
of 2 between the vortex and an antivortex enhancement 
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which would signal the Major ana state. To avoid changes 
in the background, we suggest to leave the tip at the 
same position while reversing the magnetic field. The 
field should be localized and weak enough so as not to 
overturn the chiral order parameter. 

In real three dimensional samples, zero bias peaks are 
somewhat suppressed by bulk states and surface imper- 
fections, (Ref. [22 reports a 15% enhancement above the 
high voltage background). Nevertheless, it is the ratio 
of 2 between the vortex and an antivortex enhancement 
which would signal the Major ana state. To avoid changes 
in the background, we suggest to leave the tip at the 
same position while reversing the magnetic field. The 
field should be localized and weak enough so as not to 
overturn the chiral order parameter. 

General Cases. Some difference between vortex and 
antivortex excitations is expected for any chiral sym- 
metry breaking (CSB) superconductor. The important 
questions are whether this difference is observable at 
T > Sc^ and whether it is sensitive to the existence of 
Majorana fermions. 

For a CSB superconductor with relative angular mo- 
mentum M — 1,2,... (chiral-p, chiral- (i, ...), there is 
a Majorana state in the vortex core, provided the vor- 
ticity TV obeys TV + M = 0, ±2, ±4, . . . [9 . However, in 
most cases these Majorana states vanish at the origin. 
The only exceptions are cores of antivortices which sat- 



isfy N = —M . Our factor of 2 effect will be observable 
only for this subset of cases. Notice that Px^Wy is the 
only case where the effect occurs for vorticity \N\ = 1. 

Since the core states are only sensitive to large poten- 
tial gradients, moderate disorder does not destroy the 
Majorana states. We have explicitly confirmed this ex- 
pectation numerically by solving the BdG equation 
with a white noise potential. For disorder potential fiuc- 
tuations up to order ej?, the Majorana tunneling energy 
decays with i?, and UQ{r) with r, with the same expo- 
nents as the clean system, and the peak height doubling 
signature of the Majorana states is essentially unaffected. 

Summary. We solved the BdG spectrum of Px^Wy 
vortex pair state in the spherical geometry. We show 
that even at high temperatures compared to the CdGM 
state spacing, a signature of the Majorana state remains 
when one compares the LDOS of the vortex to that of 
the antivortex. 
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